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Abstract 

The linear a model with broken U3 x U3 is compared with data on the lightest 
scalar and pseudoscalar mesons. When 5 of the 6 parameters are fixed by the 
pseudoscalar masses and decay constants one finds that, already at the tree 
level, a reasonable description for the 4 scalar masses, mixing and up to 8 tri- 
linear couplings of lightest scalars, taken as ao(980), /o(980), cj(w 500) and 
i^Q(1430). This clearly indicates that these scalars are the chiral partners of 
the vr, rj, rj' , K and strongly suggests that they like the latter are (unitarized) 
qq states. 

I. INTRODUCTION 

As is well known the naive quark models (NQM) fails badly in trying to understand 
the lightest scalars, the ao(980), /o(980), if*(1430) and the a(400 - 1200), which we shall 
here call (t(500). Therefore today most authors want to give the ao(980), /o(980) and the 
cr (500) other interpretations than being qq states. Popular alternative interpretations are 
KK bound states, 4 quark states, or for the cr, a glueball. But there is also an obvious 
reason for why the NQM fails: Chiral symmetry is absent in the NQM, but is crucial for 
the scalars. Chiral symmery is widely believed to be broken in the vacuum, and two of the 
scalars (cr and /o) have the same quantum numbers as the vacuum. Thus to understand the 
scalar nonet in the same way as we believe we understand the vectors and heavier multiplets, 
and to make a sensible comparison with experiment, one must include chiral symmetry in 



addition to flavour symmetry into the quark model. 

The simplest such chiral quark model is the linear U3xU3 sigma model with 3 flavours. 
Then we can treat both the scalar and pseudoscalar nonets simultaneously, and on the 
same footing, getting automatically small masses for the pseudoscalar octet, and symmetry 
breaking through the vacuum expectation values (VEV's) of the scalar flelds. 

As an extra bonus we have in principle a renormalizable theory, i.e. "unitarity correc- 
tions" are calculable. In fact, in the flavour symmetric {u = d = s below) limit many of the 
unitarity corrections can be considered as being already included into the mass parameters 
of the theory, once the original 4-5 parameters are replaced by the 4 physical masses for the 
singlet and octet and 0"'""'' masses, and the pseudoscalar decay constant. 

Unfortunately this over 30 years old model [0] has had very few phenomenological ap- 
plications, although important exceptions are the intensive efforts of M. Scadron and col- 
laborators 0. The reestablishment of the light and broad a has also more recently 
revitalized the interest in the linear sigma model 

II. THE LINEAR SIGMA MODEL WITH 3 FLAVOURS 

The well known linear sigma model generalized to 3 flavours with complete scalar {sa) 
and pseudoscalar {pa) nonets has at the tree- level the Lagrangian the same flavour and chiral 
symmetries as massless QCD. The U3xU3 Lagrangian with a symmetry breaking term Csb 
is 

C = ^TT[d,Ed,E^ - i/i'TrpSt] - ATrpStSSt] - A'(Tr[SSt])2 ^ _ ^-^^ 

Here S is a 3 x 3 complex matrix, = S + iP = J^^^oisa + ipa)^a/ V^, hi which 
Xa are the Gell-Mann matrices, normalized as Tr[AaA;,] = 26ab, and where for the singlet 
Ao = (2/A^/)V2i is included. Each meson in Eq. (1) has a deflnite SU3/ symmetry content, 
which in the quark model means that it has the same flavour structure as a qq meson. Thus 
the flelds Sa and pa and potential terms in Eq. (1) can be given a conventional quark line 
structure ^ (Fig. 1). 
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FIGURES 

FIG. 1. Graphical quark line representation of the A and A' terms of Eq. (1) 

The symmetry breaking terms are most simply: 

ll^SB = eaCTuu+dd + (^ss(^ss + P[det S + det S'^] , (2) 

which give the pseudoscalars mass and break the flavour and Ua{1) symmetries. The stability 
condition, that the linear terms in the fields must vanish after the shift of the scalar fields 
{Ti —>■ + V) determines the small parameters in terms of the pion and kaon masses and 
decay constants. One finds e^- = m^/yr, Gss = C^rnj^fK — m^/^)/-\/2, while (3 in the Ua{1) 
breaking term is determined by m^/ , or by + m^, . 

My fit to the scalars with the unitarized quark model (UQM) is essentially a unita- 
rization of eq.(|l|) with A ~ 16 and A' = 0, and with the main symmetry breaking generated 
by putting the pseudoscalar masses at their physical values. The model was used as an ef- 
fective theory with a symmetric smooth 3-momentum cutoff 0.54 GeV/c given by a gaussian 
form factor. Such a form factor is natural, since physical mesons are of course not pointlike, 
but have a finite size of 0.7-0.8 fm. (See the discussion in connection to Eq.(P7D below.) 
The fit included the Adler zeroes which follow from eq.(l), but only approximate crossing 
symmetry. 

Here I shall study the theory at the tree level, leaving the detailed discussion of the 
unitarization for future work. In fact, when tadpole loops are included in the unitarization 
the "unitarity corrections" to the masses should not be too large, since the tadpole loops 
partly cancel the (log A divergence in) s-channel hadron loops. One expects the corrections 
to the mass spectrum to be at most of the same order as the flavour symmetry breaking, 
since because of the renormalizability, one can in the flavour symmetric limit, include the 
unitarity corrections into the mass parameters, fixed by experiment. 

Eq. (1) without Csb is clearly invariant under S —>■ UiTiUl^ of U3xU3. After shifting 
the flavourless scalar fields by the VEV's {H —>■ H + V) to the minimum of the potential, the 
scalars aquire masses and also the pseudoscalars obtain a (small) mass because of Csb- Then 
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the A and A' terms generate trilinear spp and sss couplings, in addition to those coming 
from the Ua{^) symmetry breaking determinant term. The A term, which turns out to be 
the largest, obeys the OZI rule, while the A' and (3 terms violate this rule. 



III. TREE-LEVEL MASSES. 

It is an ideal problem for a symbolic program like Maple V to calculate the predicted 
masses, and couplings from the Lagrangian, which has 6 parameters, A, A', 13, u — d and 
s, of which the last two define the diagonal matrix V with the flavom^less meson VEV's: 
V = diag[u,d,s\. These are at the tree level related to the pion and kaon decay constants 
through u — d —< a^u^^d > /v^ = /tt/a/S (assuming isospin exact) and s —< Ggs >— 
{2 fx — /7r)/-\/2. One finds denoting the often occurring combination //^ + 4A'(ii^ + d'^ + s^) 
by and expressing the flavourless mass matrices in the ideally mixed frame: 
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+ 4:X{u^ + d^- ud) + 2(3s 



m: 



V' 



m 



m 



ss-fo 



su] 



2(3d 



eigv 



^ fl^ + 2X{u^ + d^) - 2f3s -pV2{u + d)^ 



V 



-^V2{u + d) + 4As2 j 

Jx^ + 4:\{u^ + d^ + ud) -2f3s 



^Jj? + 4A(m2 + + su) - 2(5d 



— eigv 



^ + A\'{u + d)2 + 6A(ii2 + + 2 + /3s (4A's + (5)^/2{u + d)^ 



V 



- arctan — , , ^ 
2 2A(m2 + rf2 _ 



(4A's + p)V2{u + d) 
-2V2p{u + d) 



fi' + 8A's^ + 12As 



2(3s 

2v^(4A's + /3)(M + d) 



2 4A'[(m + d)2 - 2s2] + 6A(m2 + rf2 _ 2s2) + 213s 



(3) 
(4) 

(5) 

(6) 
(7) 

(8) 

(9) 
(10) 



where eigv means the eigenvalues of the matrix which follows. 

Let us first discuss the flavour symmetric limit u — d — s. Then the pseudoscalar 
decay constants are equal fp — \f2u — f^^ = fx, while the mixing angles Qv'-stngiet _ 
^ss-f?' —54.73° and @'^-««"9'et = 0'**"'^ + 35.26° vanish, and one has 4 nondegenerate physical 
masses for the octet and singlet mesons mpsjmssjiTipojmso- Then there are also simple 



relations between the 5 model parameters A, A',/i^,/3 and u = d = s and the 4 physical 
masses and the decay constant fp: 

A = {3ml, - 2m|o - m?,s)/(12/|) (11) 

A' = (m|o - ml, - ml, + ml,)/{12fj,) (12) 

= {-3ml, + 5ml, + 4m|o)/6 (13) 

P = {ml,- ml,) /{3V2 fp) (14) 

u = d = s = <auu >=< cr^d >=< cTss >= fp/ V2 (15) 



It is obvious that we can reparametrize the theory in terms of these physical quantities, 
which can be kept fixed in the renormalization, as long as fiavour symmetry is exact. (If one 
choses the same tree-level values for the parameters A,/i^,/? as found below in Eq.(|l6D and 
A' = 1 but foTU = d = s the average value, 75.05 MeV, or fp = 106 MeV, one would have 
mps = 384 MeV, mpo = 956 MeV, mss = 1086 MeV, mso = T41 MeV.) Thus these masses 
can be thought of as already "unitarized" . On the other hand the original parameters and 
induced tri-meson couplings will be renormalized from the tree level values. 

Now breaking the fiavour symmetry {s u = d) we have only one more parameter, 
given by fx — fn, and it is evident that this breaking splits the degeneracy in the mass 
spectrum from 4 independent masses to 8 masses and generates two mixing angles. Thus 
one gets several tree- level predictions in particular for the scalars (Table I). Of course now we 
expect these tree-level predictions to receive corrections from the unitarization, but for small 
symmetry breaking (experimentally (s — n)/s ~ 0.308) one would expect these corrections 
not to be larger than this, i.e. <30%. 

We can fix 5 of the 6 parameters, leaving A' free, by the 5 experimental quantities from 
the pseudoscalar sector alone: m^, m^, + m^/, f^, = 92.42 MeV and fx = 113 MeV 0, 
which all are accurately known from experiment. One finds that at the tree level 

A = 11.57, if = 0.1424 GeV^ P = -1701 MeV, u = d = 65.35 MeV, s = 94.45 MeV. 

(16) 

The remaining A' parameter changes only the a and /o masses and their trilinear couplings, 
not those of the pseudoscalars. Gavin et. al. P] calls for a group theoretical reason for this 
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simplificati on. In fact, graphically it is easy to see that the A' term can only break the OZI 
rule for the scalars, which can couple to the vaccuum, but for the pseudoscalars the A' term 
must leave the OZI rule intact. Therefore A' and /z affects the pseudoscalars only through 
the combination /i^ = /x^ + 4A'(-u^ + rf^ + s^). 

It is of some interest that the simple Gell-Mann-Okubo mass formula for the mixing does 
not give the same result for the mixing angle between t] and rj' as our model. This is because 
there is flavour symmetry breaking also in the anomaly terms f3d and f3s in Eqs. (3-10) above. 
E.g. for the octet pseudoscalar mass one gets from Eqs. (3-10) (4m|^ — m^)/3 = 542.5 MeV 
and a mixing angle of B**~^' = —12.7°. This would be closer to the conventional mixing 
angle (—10° to —23°) than our model —5.0°, where the octet rjg mass is 566.1 MeV. 

Some of the couplings of a and /o depend sensitively on A', since A' changes the small 
ideal mixing angle, (f)'^'^~^° . It turns out below that A' must be small, compared to A, in order 
to fit the tri-linear couplings. By putting A' = 1 one gets a reasonable compromise for most 
of these couplings. With A' ~ 3.75 one almost cancels the OZI rule breaking coming from 
the determinant term, and the scalar mixing becomes near ideal (for A' = — /3/(4s) = 4.5 
the cancellation is exact). 

As can be seen from Table I the predictions are not far from the experimental masses 
taken as ao(980), /o(980), i^'Q(1430), and cr(500). For a discussion of the existence of the 
cr(500) see my recent Frascati talk which also includes some preliminary results of the 
present paper. Considering that one expects from our previous UQM analysis of the scalars 
and our discussion above that unitarity corrections can be up to 30% , and should go in 
the right direction compared to experiment, one must conclude that these results are even 
better than expected. Similar mass analyses as in Table I have been done in Refs. 
although with somewhat different input data. 

IV. TRI-LINEAR COUPLINGS AT THE TREE LEVEL. 

The trilinear coupling constants follow from the Lagrangian after one has made the shift 
^ H + V . The most important spp couplings are at the tree level, when expressed in 
terms of the original parameters: 
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= 4A(d -u + s)-2p 






(17) 
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= sin 0^""--'^° a/2 [A (4m - 2s) + 4A'(m + d) + (3] - A cos 0'^" 


■^"[A(m 


- 2s) - 2A's] 


(23) 




= cos0''"-'''2V2A(m + d)- 2(3 sin (j)''-"' 






(24) 




= sin (f)''-'''2V2X{u + d) + 2p cos (f)''-''' 






(25) 


QaoK+K- 


= V2[A(4m - 2s) - 4A'(d - u) - p] 






(26) 



In fact these can be written in more useful forms in terms of the predicted physical masses 
and mixing angles and decay constants: 
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(36) 



In table II 8 different spp couplings are compared with quoted experimental numbers. In 
some of the channels of table II the resonance is below threshold and the widths therefore 
vanish at the resonance mass. However, the coupling constants have recently been deter- 
mined through a loop diagram from — > KK — > 77r7r and — > KK —>■ ■jttt) (albeit in a 
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somewhat model dependent way) by the Novosibirsk group 0,|T^]. For channels where the 
phase space is large, it is important that one includes a form factor related to the finite size 
of physical mesons. In the ^Pq quark pair creation model a radius of ~0.8 fm leads to a 
gaussian form factor, as in the formula below, where ko ^ 0.56 GeV/c (as was found in the 
UQM 0). Thus the widths are computed from the formula: 

As can be seen from table II most of the couplings are not far from experiment. Only 
the /o TTTT and ao —>■ ni] couplings and widths come out a bit large, but these are very 
sensitive to higher order loop corrections due to the KK threshold, and /o —>■ tttt is extremely 
sensitive to the scalar near-ideal mixing angle and A'. If one choses A' = 3.75 this mixing 
angle nearly vanishes {(f)'^^~f° = —3.0°) together with the /o tttt coupling (c.f. eq. (|32D ). 
From our experience with the UQM |^ the ao KK peak width, when unitarized, is 
reduced, because of the KK theshold, by up to a factor 5. Therefore one cannot expect 
that the tree level couplings should agree better with data than what those of Table II do. 
After all, this is a very strong coupling model (A = 11.57, leading to large gf /Atc) and higher 
order effects should be important. 



V. CONCLUSIONS. 

In summary, I find that the linear sigma model with three flavours, at the tree level, 
works much better than what is generally believed. When the 6 model parameters are fixed 
mainly by the pseudoscalar masses and decay constants, one predicts the 4 scalar masses 
and mixing angle to be near those of the experimentally observed nonet ao(980), /o(980), 
(t(500), -^"0(1430). Also 8 couplings/widths of the scalars to two pseudoscalars are predicted 
reasonably close to their presently known, rather uncertain experimental values. The agree- 
ment is good enough considering that some of these are expected to have large higher order 
corrections. The model works, in my opinion, just as well as the naive quark model works 
for the heavier nonets. A more detailed data comparison would become meaningful, after 
one has included higher order effects, i.e. after one has unitarized the model, e.g., along the 
lines of the UQM g. 
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Those working on chiral perturbation theory and nonhnear sigma models usually point 
out that the linear model does not predict all low energy constants correctly. However, one 
should remember that the energy regions of validity are different for the two approaches. 
Chiral perturbation theory usually breaks down when one approaches the first scalar res- 
onance. The linear sigma model, on the other hand, includes the scalars from the start 
and can be a reasonable interpolating model in the intermediate energy region near 1 GeV, 
where QCD is too difficult to solve. 

These results strongly favour the interpretation that the ao(980), /o(980), cr(500), 
/rQ(1430) belong to the same nonet, and that they are the chiral partners of the tt, r}, 
rj', K. If the latter are beheved to be unitarized qq states, so are the hght scalars ao(980), 
/o(980), (7(500), Kq{1A30), and the broad (7(500) should be interpreted as an existing reso- 
nance. The cr is a very important hadron indeed, as is evident in the sigma model, because 
this is the boson which gives the constituent quarks most of their mass and thereby it gives 
also the light hadrons most of their mass. Therefore it is natural to consider the cr(500) as 
the Higgs boson of strong interactions. 

Acknowledgement. This work is partially supported by the EEC-TMR program Contract 
N.CT98-0169. 
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TABLES 

TABLE L Predicted masses in Me V and mixing angles for two values of the X' parameter. The 
asterix means that m-,^,mK and rn^ + m^, are fixed by experiment together with f^^ =92.42 MeV 
and fx =113 MeV. 



Quantity 


Model A' = 1 


Model A' = 3.75 


Experiment 


"T-TT 


137*) 


137*) 


137 i 


TTlK 


495*) 


495*) 


495 [| 




538*) 


538*) 


547.3 i 




963*) 


963*) 


957.8 i 


QTj' —singlet 


-5.0° 


-5.0° 


(-16.5±6.5)° § 


ruao 


1028 


1028 


983 i 




1123 


1123 


1430 II 




651 


619 


400-1200 @ 




1229 


1188 


980 [| 


Qcr— singlet 


21.9° 


32.3° 


(28-18.5)° i 
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TABLE II. Predicted couplings GeV^) , when X' = 1, compared with experiment 

and predicted widths with experiment (in MeV). (We have used isospin invariance to get the sum 
over charge channels, when there is data for one channel only.) The predicted /o vrvr width is 
extremely sensitive to the value of X' (for X' = 3.75 it nearly vanishes) and unitarity effects as 
discussed in the text. Also the oovrr/ coupling is very sensitive to loop corrections due to the KK 
threshold. 



Process 


V 3? 

Z^i 47r 

in model 


V 3^ 

in experiment 


model 


experiment 


a vr+vr^ + 7r°7r° 
a K+K- + K^K^ 
fo vr+vr- + 7r°7r0 
/o ^ K+K- + K^K° 
at TT+r? 
a+ ^ K+K'> 


7.22 
0.28 
2.17 
0.16 
1.67 
6.54 
2.29 
2.05 


«o 

1.95 [11| 
0.004 111] 

0.7651°:?° i 

4.26^Pi2 |12 

0.57 [12| 
1.34l°:i |12 


1 


678 

13 
574 



see text 


273 see text 



278 ± 23 
< 26 [10| 
300-1000 II 


40 - 100 [| 


50 - 100 [| 
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